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Abstract 

We study the effect of a sharply localized magnetic field on the mobility of electrons in a strip 
(ribbon) of graphene sheet. The magnetic field is taken as a single and double delta type localized 
functions, which are the zero width limit of gaussian fields considered at a later stage. For both 
field configurations, we investigate the scattering behavior of electrons and evaluate analytically 
and numerically their transmission and reflection coefficients. The possibility of spacial confinement 
due to the inhomogeneous field configuration is also investigated. 
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1 Introduction 



Graphene, a single sheet of pure carbon atoms, has been the subject of intense focus and research for 
the last few years due to its unique and amazing properties that could make it more suitable for future 
electronic devices and possibly replace silicon-based semiconductors in electronics [TJ [2]. Graphene 
nanoribbons, in particular, are of interest because they exhibit a conduction bandgap that arises due 
to size confinement and edge effects [3j-[llJ. From the electronic applications point of view, graphene- 
based materials with gap are very desirable. This lead to substantial theoretical and experimental 
efforts to generate and control the energy gap in graphene based devices. One approach is based on 
quantum confinement as in quantum dots and nanoribbons, where it was shown that the energy gap 
value increases with decreasing nanoribbon width |12j . 

Theoretical studies have suggested that graphene nanoribbons could have interesting magneto- 
electronic properties, with predicted very large magnetoresistance [HdH]. Experimental observation 
of a significant enhancement in the conductance of a graphene nanoribbon field-effect transistor by a 
perpendicular magnetic field has been reported recently [14]. They have observed an extraordinarily 
large tunable magnetoresistance in graphene nanoribbon-FET devices and an increase in conductance 
by a factor of over 10,000 was demonstrated at 1.6K. This magnetoresistance can be readily tuned 
by the gate voltage and source-drain bias, with enhancement reaching maxima near the edges of the 
conduction bandgap. These experimental findings clearly demonstrate that the graphene nanoribbons 
exhibit interesting magneto-transport properties, with the possibility of opening up exciting opportu- 
nities in magnetic sensing and a new generation of magneto-electronic devices. 

On the other hand, more efforts have been employed to understand the scattering behavior of 
fermions in graphene. In particular we cite [15] where the transmission through magnetic barriers in 
graphene-based nanostructures. Several particular cases were considered: a magnetic step, single and 
double barriers, and ^-function barriers. A separate class of magnetic-barrier structures are those with 
inhomogeneous magnetic- field profiles, such that the average magnetic field vanishes, which can be 
realized by nanostructured ferromagnetic stripes placed on top of the graphene layer. Quantum bound 
states that are localized near or in the barrier are predicted for a magnetic step and some structures 
with finite-width barriers but none for 5-function barriers. When a bound state is localized close to 
the barrier edge, it has a nonzero velocity parallel to this edge. The transmission depends strongly 
on the direction of the incident electron or hole wave vector and gives the possibility to construct 
a direction dependent wave vector filter. In general, the resonant structure of the transmission is 
significantly more pronounced for (Dirac) electrons with linear spectrum than for the usual electrons 
with a parabolic spectrum. 

Motivated by different results on the scattering behavior, we study fermions in graphene nanorib- 
bon subject to a sharply localized magnetic field. Indeed, we show that the solution space of the 
present system is the union of the two solution subspaces governed by two different kinetic balance 
equations. These two solution subspaces are not identical and their union constitutes the whole space 
of solution of our problem. To resolve the issue of which energy domain is selected by each kinetic 
balance relation, we perform a unitary transformation on the fundamental Dirac equation. Such a 
transformation does not affect the energy spectrum or the physics of the problem. We choose a rota- 
tion by 7r/4 about the y-axis, U = e*4 <72 . After matching the eigenspinors at the boundary conditions, 
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we evaluate the reflection and transmission amplitudes in terms of the energy and magnetic field. A 
detailed study of their coefficients is presented as well. 

Subsequently, we consider a magnetic field normal to the graphene sheet of a bell shape, i.e. 
B(x) = zBosech 2 (iJ,x) where -Bo is the strength of the field and stands for its range. We show 
that the solution space of the corresponding Hamiltonian splits into three regions depending on the sign 
of fci?o and the energy range, with k wave vector along the y-axis. Using the corresponding current 
density, we evaluate the reflection and transmission coefficients, which reduce to those obtained in 
section 3 for \x = 1000. Finally, we treat the case of a double step pseudoscalar potential and also 
double Gaussian type magnetic field. 

Our manuscript is organized as follows. We present our theoretical model by setting the general 
form of potential in section 2. As a first illustration of potential type, we consider 5-type magnetic 
field in section 3. After getting the energy spectrum solutions, we study the scattering problem by 
evaluating explicitly the reflection and transmission coefficients and underline their physical properties. 
Because of 5- type magnetic field is a limit case of a general form, we study Dirac fermions in a Gaussian 
magnetic field, in section 4. These studies suggested to deal with an integrating case that is Dirac 
fermions in a double 5-type magnetic field in section 5 and the corresponding double Gaussian type 
magnetic field is considered in appendix. We conclude our work in the final section. 



2 Theoretical model 

In graphene, the low energy linear electronic band dispersion near the Dirac points gave rise to charge 
carriers (electrons or holes) that propagate as if they were massless fermions with speeds of the order 
of 10 6 m/s rather than the speed of light 3 x 10 8 m/s. Hence, charge carriers in this system should be 
described by the massless Dirac equation rather than the usual Schrodinger equation. However, we 
should stress that such a continuum approximation remains valid up to wave vectors of the order of 
the inverse lattice spacing. Under this approximation we can describe our charge carriers by a gauge 
invariant Dirac Hamiltonian in 2 + 1 space-time dimensions with minimal coupling to the three- vector 
potential = (Aq,A) = (V, A x ,A y ) as follows 

H = mj + V + j j ■ (iV + A) (1) 

where 70 = 03, 7 = {^ x ,ly) = i{&i,&2) and {ai]\ =l are the 2x2 Pauli matrices. Adopting the 
conventional relativistic units H = c = 1, we obtain [16] 



As we claimed before, for a system made of graphene the charge carriers are massless m = and 
moving with the Fermi velocity (about 10 6 m/s). Moreover, we specialize to a physical configuration 
corresponding to a single graphene sheet in a magnetic field with vector potential A = yW(x). Thus, 
the associated stationary Dirac equation becomes 



0. (3) 
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We can separate variables and write 



(4) 



where x{u) = e lky anci ^ is a real parameter that stands for the wave number of the excitations along 
the y-axis. The resulting reduced equation is given by 



dx 



+ W + k 



dx 



+ W + k 



\<P- J 



(5) 



In the present representation W plays the role of a pseudoscalar potential coupling. The above 
equation is equivalent to the following set of equations 



±^ + W + k 
dx 



± 



ax z ax 



± 



0. 



(6) 
(7) 



The structure of the Schrodinger-like equation ^ has the features associated with supersymmetric 
quantum mechanics (SSQM) |17j . The two potential partners in SSQM are V 2 =t V', where V is the 
superpotential and in our case we have V = W + k. The relativistic extension of SSQM exhibited 
in some models that are governed by the Dirac theory is due to the fact that the associated Dirac 
Hamiltonian is an element of the 4-dimensional superalgebra £7(1/1), which is a Z2 graded extension 
of 50(2,1) Lie algebra [18J. It is interesting to note that the corresponding massive (1+1) Dirac 



equation 



m — e 



dx 

—m 



+ W 

£ 



/ ^ 



(8) 



results in the following Schrodinger-like equation 



d dW 2 

- —j + W =F ~, \~ m 

ax z ax 



r 







(9) 



which has the same structure associated with SSQM as in with the mapping W — > W+k. Therefore, 
the solution of the two problems are similar with k playing the role of the mass m. We have encountered 
this phenomenon before, where space reduction (due to compactification or symmetry) transforms the 
2D massless problem into an effective ID massive problem [19] . Note that the complete solution space 
of the Dirac equation ^ is the union of the two solution subspaces of equations ^ and Q with 
the top and bottom signs, respectively. Below, we identify the energy regions associated with each 
subspace. 



3 Dirac Fermions in a c^-type magnetic field 

Here we consider the zero range and infinite strength limit where the magnetic field takes the form 
B{x) = zBq8(x/o), with a is a length scale parameter. The main advantage of this configuration is 
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that it is amenable to analytical treatment and interesting results. Thus, the single-step pseudoscalar 
potential is given by 

W(x) = aB @(x/a) (10) 



where G(x) being the Heaviside step function, such as 



Q(x) 



1. 
0. 



x > 
x < 



mi 



Figure 1 shows the schematics of a laboratory set-up where a coil of a rectangular cross section with 
very narrow width will result in the necessary current sheet that should produce a magnetic field 
approximating that of our problem. The graphene ribbon is the strip that penetrates the coil along 
the x-axis in the figure: 




Fig. 1: Schematic of a set-up where a coil of a rectangular cross section with very narrow width produces 
a magnetic field approximating that of our problem. 

To study the scattering behavior in the present problem, we first determine the plane wave solutions 
associated with both regions, dbsc > 0. The two components of the wavefunction satisfying equations 
^ and with the top sign are given by 

1 / d 



e V dx 



+ u)±\<p 



dx 2 



2 2 

+ u4 - e 



if + = 



(12) 
(13) 



with uj + = oBq + k and w_ = k. It is easy to show that the solution of these equations reads 



<p(x) = A 



K±X 



+ B 



-K±X 



Ci± 



01 



where the parameters read as 

k± = Ju\-£ 2 , a± = (w± -k±)/e, P± = (u± + k±) /e. 



(14) 



(15) 
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On the other hand, considering the lower sign in equations ([6]) and ([7]) for ±x > we obtain 

, 1 / d \ _ 

e \ ax J 
d 2 1 

Uj — 

+ a4 - er (fi =0 



dx 2 



(16) 
(17) 



which can be solved to give 



<p(x) 




K±X 



+ B 




-K±X 



(18) 



In summary, the oscillatory solutions (for both cases and for ±x > 0) are obtained for all |e| > \oj±\ 
while evanescent wave solutions will result if |e| < \jj±\. Thus the quantity \ui±\ plays the role of an 
effective mass with different asymptotic values at ±x > 0. That is, our massless problem is equivalent 
to a space dependent massive problem whose effective mass depends on both the sign and strength of 
the localized magnetic field and the value of component of the wave number along the y-direction. 

We need to reassert that the total solution space of the problem is the union of the two solution 
subspaces ^ and ([7]). These two solution subspaces are not identical and their union constitutes the 
whole space of solution of our problem. To resolve the issue of which energy domain is selected by 



which kinetic balance relation (12) or (16), we perform a unitary transformation on the fundamental 



Dirac equation (pi). Such a transformation does not affect the energy spectrum or the physics of 

— ■ TL 

the problem. We choose a rotation by 7r/4 about the y-axis, U = e . Thus, the transformed 
Hamiltonian and wavefunction read 



H 
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UHU ] 
Utp 



4(72 dx + + k ^ 3 



(19) 



then ([5]) becomes 



H-el\(p = 

leading to the following new kinetic balance relation 



d 



1(72 dx + ( W + k } (X3 



el 



d 



e±(W + k)dx 



which is not valid for 



e = =F (aB Q + k) = =Fw^ 



(20) 



(21) 



(22) 



Thus, in the case when aB^ is a positive constant, this relation with the top/bottom sign corresponds 
to the positive/negative energy subspace^J Now, since equation (20) results from a unitary transfor- 



mations of (12) and (16), then the same applies to those equations as well. That is, equations (12) 



and (16) are valid for positive and negative energy, respectively. Figure 2 shows the potential con- 



figuration of the problem. Dark/light grey areas correspond to negative/positive energy continuum 
where the solution is oscillatory of the form e ±lKX . The white areas correspond to bound states where 



2 If k is taken negative in equation (4), then equation (14) with the top/bottom sign corresponds to negative/positive 
energy subspace if W+k becomes negative. 
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the solutions are exponentials, e ±KX . In the subsequent developments, we will only be considering 
scattering solutions with a normalized incident beam from left. Thus, the relativistic energy must be 
e > \k\ and the solution space is built only from the positive energy solutions (14). 



1 => 



l+|<u + | 
■ aB„ 




Fig. 2: The potential configuration of the problem. Dark/light grey areas correspond to negative/positive 
energy continuum. The white areas correspond to bound states where the solutions are exponentials. 



The total solution space splits into two regions depending on the energy range as follows. For 
k < e < u>+, scattering with total reflection will take place. It should be noted that there is no Klein 
tunneling in this problem due to the absence of vector interaction. With a normalized incident flux 
from left, the solution of equation ^ is obtained from (14) as 



A Kl 



-K+X 



x < 
x > 



(23) 



where k = \n—\ 
reflection amplitude 



V e 2 — k 2 and a = (k — %K)je. Continuity of the wavefunction at x = gives the 
(3 + — a clBq + k + + in 



R=- r =~ Z ( 24 ) 
p + — a* clBq + k + — ik 

and the normalization parameter A = 1 + R. It is clear that we have complete reflection in this case 

since \R\ = 1. It is also worth mentioning that the probability density in the positive x-direction is 

given by 

p{x) = \A\ 2 e- 2 \ R +\ x {l + \P + \ 2 ) (25) 



which shows that the penetration depth, A, of the Fermion into the barrier region defined by 

A 1 



2|At+| 2 2y/(aB + k) 2 - e 2 



(26) 



decreases very rapidly with increase in the magnetic field strength. This result indicates the possibility 
to confine the fermions in any region between two <5-like magnetic fields generated by a symmetric 
double step pseudoscalar potential, this situation will be treated in section 5. 
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Now, let us move on to analyze the scattering situation when e > uj + so that partial reflection will 
take place. Again, with a normalized incident wave from left, the solution of ^ is obtained from (14) 

as 

, . / Q)e^+RQ,)e-^, x<0 
^ X) = { T *>0 (27) 



where 7 = \k+\ = ye 2 — uo\ and j3 = — £7) /e. Again matching the spinor wavefunctions at x = 
gives the reflection and transmission amplitudes 

a = -4=4 r=^. (28) 

One can verify the flux conservation relation 

\R\ 2 + | 7 /k| |T| 2 = 1. (29) 

Figure 3 is a plot of the reflection and transmission coefficients as a function of the energy for several 
values of the magnetic field strength Bq. 




It is also worthwhile to study the transmission coefficient as a function of the magnetic field strength 
Bq. For this purpose we can write explicitly |T|, in this regime, as follows 



k 2 



\ (aS ) 2 + (V^¥ + ^e 2 - (aB + k) 2 



(30) 



which shows clearly that the transmission is unity at -Bo = and reduces to zero for strong field. 
Figure 4 shows the behavior of the transmission and reflection coefficients as a function of the applied 
magnetic field. We would like to note that the transmission vanishes for energies below the critical 
energy value of e* = Bq + k, as shown in figure 3. A similar remark can be made regarding figure 4. 
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Fig. 4: Reflection and transmission amplitudes as a function of strength of the field for e = 3,4,8 in units 
of a~ 2 (left to right curves) and k = 2. 



4 Dirac fermions in a Gaussian magnetic field 

We consider a magnetic field normal to the graphene sheet of a bell shape. Specifically, we take 
B(x) = zBq sech 2 (/xx) (Fig. 5) |20l[2T], where Bq is the strength of the field and // _1 stands for its 
range. The corresponding vector potential (Fig. 5) is 

A{x) = yW{x) = y (B /m) tanh(/xx) (31) 

which turns equation Q into 

/ n_ \ 2 1 

ip ± = 0. (32) 



d 2 „ ^ , o x f B, 



" 2 ^B) (1 — tanh 2 /tx) + ( — - tanh/xx + A;J — e 2 
Note that as x — > ±oo, this equation reduces to a simple form 

£+4-*)^-° (33 ) 

where wj- = k ± Bo///. This equation is similar to that of a massive ID problem with asymptotic 
masses m = |cj-i-| at x — > ±oo. The solution space of (32) splits into three regions depending on the 
sign of kBo and the energy range. If fcBo > 0, then |u>_| < and these energy regions can be 
defined as follows 

• |e| < \u}-\: Bound states with the boundary condition lim (p(x) = 0. 

x— >±oo 

• \uj-\ < \e\ < Scattering with total reflection to the negative x-axis (note that there is 
no Klein tunneling in this problem due to the absence of vector interaction). The boundary 
conditions in this case are: lim ip(x) = 0, lim (p(x) = e 1K ~ x ( * ) + lZe~ 1K ~ x I i ), where 



k± = e 2 - lo±, p± = I (uj± — ik±) and the reflection amplitude is such that \1Z\ = 1. 
|e| > \u)+\: Scattering solutions occur with the boundary conditions: lim (f(x) = Te 1K+x ( 1 



lim tp(x) = e 1K ~ x (p_) + Tie. 1K ~ X ( i ), where is T the transmission amplitude such that 
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\1Z\ 2 + \p + / p-\ \T\ 2 = 1. However, if kBo < 0, then \oj-\ > \uj + \ and the above subspaces should 
be interchanged with u + f-> w_ and x — > —x. 




WW 
0.002*— 



^=1000 




Fig. 5: Hyperbolic tangent pseudoscalar potential profile along with the associated magnetic field for the 
strength B = 2. 



To be more specific, let us consider the scattering problem in the negative x regions, x < 0. Making 



the change of variable z 



equation (32) becomes 



4&V 



d ( d 



dz V dz 



TB 



42 



{l-zf 



B 1 + z i ■■> 

■ — Vk) +e z 

1 — z 



ip^ = 0. 



(34) 



Now, we solve (34) for and then deduce <p using Q. The general solution of (34) can be expressed 
in terms of the hypergeometric functions [22j as follows 



where we have set h! 



Az h '(l - z) x F (ti - v + \U + v + A, 1 + 2h; z) 
Bz- h '{\ - z) x F {-h! - v + A, -ti + v + A, 1 - 2ti; z) 



(35) 



2k ' 



■a, and \ = \ 



To study the scattering 

problem we need to investigate the asymptotic behavior of (35) when x — > — oo, which corresponds to 



/fe2+4B +4(^l) 2 
2k ' 



z — > 0. Using the asymptotic behavior of the hypergeometric functions we obtain for ip + 

lim ip + = Ae ik - X + Be~ lk - X (36) 
with A;_ = 2ifj,\h'\. From ^ the other component, ip~ , is given by 



e \dx /i 



Substituting equation (36) in (37) gives the asymptotic behavior of ip 

-ik—X 



lim <p = Arj + e lk - x + Brj-e~ 

x—>—oo 



(37) 



(38) 



-i u)+±2uh' 

where r/± = — — t — . 



Next, we consider the solution of equation (32) for x > 0. The analysis of the solution can be 



simplified by making the substitution z 



-2[ix 



. Equation (32) becomes 



4k 2 z- 



d d 



dz V dz 



4z 



l-zf 



Bo 1 + z \ 

— \-k] +e z 

p 1 — z 



ip ± = 0. 



(39) 
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The general solution to the above is 

<p+ = Cz h (l-z) x F(h-v + \,h + v + \, l + 2h;z) 

+ Dz~ h (l - z) x F (-h - v + A, -h + v + A, 1 - 2h; z) 



(40) 



where h 



and A 



l 



^+4B +4(ga) 2 



2k ' " 2k " iJJa 2 2fc • 

Also as x — > — oo, z — > 0. Therefore in order to have a plane wave traveling to the right as x — > oo, 
we require that \e\ > |w+| and C = so that as x — > oo we have 



Dz~ h {l - z) A F (—h -v + A,—h + v + \,l — 2h; z) 



and the limit 



lim (p 

X— > + oo 



ik+x 



with = 2i/x|/i|. The other component is 

lim ip~ 



x—t+oo 



D(3e 



ik+hx 



(41) 
(42) 

(43) 



where f3 = 2 ^ tfe + aj + _ 

For one-dimensional scattering problems, the particles in the beam are in plane wave states with 
definite momentum. Given the wave functions relevant to incident, reflected and transmitted beams, 
one may calculate the corresponding current densities [23J . The reflection and transmission coefficients 
are defined to be the ratio of the reflected current density j re to the incident current density j m and 
transmitted current density j tr to the incident current density, respectively, leading to 



R 



\f 



V-\ i ,2 

r \ f\ 



\J\ \v+-v+\ 

where the conservation law R + T = 1 is obeyed. 



J' 



1/3-/8*1 | ,2 
i r l r l 

\v+-v+\ 



(44) 



1 
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* 
1 
1 
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T — 
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(b) 



T 
R- 



Fig. 6: Transmission as a function of energy with k = 1, Bo = 2 and (a) n = 3, (b) \x = 1000. 

We notice that for the Hyperbolic tangent magnetic field (figure 5), figure 6 show clearly that for the 
value of [i = 1000 we obtain the same result for the transmission and reflection coefficients as those 
generated in section 3 for a single step pseudoscalar potential. 

We can also characterize the transmission amplitude by giving the corresponding phase shifts 5. 
Indeed, in the case of < |e| < \oj+\ one obtains 



t = \t\ e 



id 



(45) 
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The calculations lead to the scattering phase shift in the form 

h 



arctan 



h' 



h" 



e 2 -u?_ 



2k 



(46) 



A complementary information about the phase shift can be obtained by studying the trajectories in the 
complex transmission amplitude plane where full reflections will appear whenever these trajectories 
pass through the origin. 

5 Dirac fermions in a double 5-type magnetic field 

Motivated by the above analysis, let us investigate the scattering of Dirac fermions in a double 5- type 
magnetic field. This can be described by the double-step pseudoscalar potential 



W{x) = Bi[l- Q(x + a)] + B 2 Q(x - a) 



(47) 



presented in figure 7 with a defined to be positive real number. This pseudoscalar double step poten- 
tial generates a magnetic field B = B\5{x + a) + B 2 5(x — a), B\ and B 2 are the strengths of the <5-type 
magnetic fields at x = —a and x = +a, respectively. 



Fig. 7: Double-step pseudoscalar potential. 



Using (47), we show that the general solutions of the Dirac equation ([7]) for the case of 
can be written as 



£ > W4 




IK — X I 1 



A+e iK ~ x 



+ B_e- 1K - 2 
+ B+e- iK - 3 



to. 



1 I ' 



c+e ik 2 x h+\ +D+e 



i k2X II 



E_e iK + x i + + F_e-' 1K + X (l« 



E,e 



i K+ x (P 



e < oj + 
e > —u + 

e < UJ + 

e > —oj + 
e < oj + 
e > —uj + 



(48) 



where 7+ = (ik2 — k)/e, n± = w e 2 — u±, k 2 = \/e 2 — k 2 , u + = B 2 + k and o;_ = B\ + k. These 
wavefunctions will be used to study the scattering and the bound states for the present case. 

The continuity of the wavefunction at x = ±a gives the reflection amplitudes = Bzp/Azp and the 
transmission amplitudes t T = E T /A T for the negative (top sign) and positive (bottom sign) energy 
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solutions. More precisely, we find 

a -ia(-2k 2 +K^+K + ) u 

+ = K - K2 f49) 

+ e 4iak 2(K„-k2TiB 1 )(k2-K+=FiB 2 ) + (K- + k2TiBi){k2 + K + ±iB 2 ) y ' 

e~ 2iaK - [e 4iak *(K_ +k 2 T iB 1 )(k 2 - k+ T iB 2 ) + («- -k 2 ± igiKfcg + k+ ± iB 2 )] 

Tt e 4iafe 2 ( K _ _ k2 ^ iB{){k2 ~ K + -F iB 2 ) + (k- + k 2 T iBi){k 2 + k+± iB 2 ) ' ( ' 



The result of the passage to the limits, B\ — > (k 2 —)■«;_) and a — > 0, in (49) and (50) coincides with 
the transmission and reflection amplitudes in ( |28| ) of the Dirac particle scattered by a pseudoscalar 
single step potential. If we take the limit a — > in the above equations we get, by replacing u + and 
uj- by uj + = B 2 — k and w_ = B\ — k, respectively, and Bq by B\ — B2, the transmission and reflection 
amplitude for the Dirac particle scattered by a double-step pseudoscalar potential W(x) M7| . For 
the special case .E>2 = — = 0.002 the transmission and reflection coefficients coincide with those 
obtained in the case of scattering of Dirac fermions by a single step pseudoscalar potential plotted in 
figure 6 (b). 

Given the wavefunctions associated with the incident, reflected and transmitted beams, one may 
calculate the corresponding current densities according to [23]. Figure 8 illustrate the configuration 
of the solution space of our problem. Dark/light grey areas correspond to negative/positive energy 
continuum where we have evanescent or transmitted waves (EWs/TWs). The white areas correspond 
to the bound states (BS). The reflection and transmission coefficients for the positive and negative 
energy solutions are 



which obey the conservation relation 



, &(«+) U |2 1 fK0 \ 



Figures 9 and 10 show more clearly that, for the confinement by a double-step pseudoscalar potential, 
the EWs, TWs and BS always appear for the energy that < \e\ < \u+\ , |e| > \uj + \ and |e| < \uJ-\, 
respectively. In figures 9 (a), (b) and 10 (b) the transmission resonances appear for B\ and B2 7^ 0. 
From figure 10 (a) one readily notice that, for the value of B2 = 0, the double-step pseudoscalar 
potential reduces to a one step with W{x) = B\(l — 0(x)). We can easily find the EWs, TWs and BS 
for this later representation by replacement the u + and cj_ by w + = k and = B\ — k respectively. 
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Fig. 8: The potential configuration associated with the double step pseudoscalar potential. Dark/light 
grey areas correspond to negative/positive energy continuum. The white areas correspond to bound states 
where the solutions are exponentials. 




Fig. 9: Reflection and transmission coefficients as a function of energy for B2 = —B\ = 0.5,2.0,3.0,5.0 
(left to right curves) with k = 1. (a) for a = 1 and (b) for a = 2. 




Fig. 10: The transmission coefficient as a function of energy for a = 2 and k = 1. (a) for the Dirac particle 
scattered by a pseudo-scalar one step potential when B2 = and B\ = —0.5,-2.0,-3.0,-5.0 (left to 
right curves), (b) for our the case of B2 = B\ = —0.5, —2.0, —3.0, —5.0. 



We wish to study the bound states of a double-step pseudoscalar potential for the case |e| < 
and for B\ = B2. From equation (48) we must have = F_ = in order to obtain normalizable 



wavefunction in the region |:r| > a/2. We now have the full solution, given our assumption of particles 
incident from the left 



93 



B_e 



K—X I 1 

a. 



E_e~ K + x (i 



-i k2X 1 1 



x < —a, 
\x\ < a, 
x > a. 



(53) 



This means that the solutions separate into even parity and odd parity states. We could have guessed 
this from the potential. The even states (C_ = = EL) have the (quantization) constraint on 

the energy that 

ak 



1 



ak 

— v + v + 
e e 



tan(e) 



(54) 
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and the odd states (C_ 



-E-) have the constraint 



ak ak 
— v + v H 

e e 



cot(e) 



(55) 



where we have used v 



Ba 
2e ■ 



These are transcendental equations, so we will solve them graphically. 



The plot in figure 11 shows the left hand side of the transcendental equation along with the tangent 
for even parity states and minus cotangent for odd parity states, respectively. The intersection of 
these curves gives the bound states. From this figure we see that there are two bound states for even 
parity solution and two for odd parity solutions. The wider and deeper the step, the more bound 
states become available. 




Tan(e) 

Cot(e) 



Fig. 11: The double-step pseudoscalar potential bound states for k = 1 and a = 1. The solid line represents 
tan(e) and the dashed line represents — cot(e). The dotted lines represent the left-hand side of equations 
d54h and (pi). 



6 Conclusion 

We considered Dirac fermions in graphene nanoribbon subject to a sharply localized magnetic field. 
After setting the theoretical model we obtained the energy spectrum solutions for a sharply local- 
ized 5-type magnetic field. After matching the wavefunction at the point x = 0, the reflection and 
transmission amplitudes are derived in the standard way. A study of the transmission coefficient as a 
function of the strength of field -Bo was done, the transmission is unity for Bq = and reduces sharply 
for strong field. 

The general case of a bell shape magnetic field normal to the graphene sheet, namely B{x) = 
zBosech 2 ([ix) was studied, Bq being the strength of the field and stands for its range. It was 
shown that the solution space of the corresponding Hamiltonian splits into three regions depending on 
the sign of kBo and the energy range, k is the wave vector along the y-axis. To obtain the reflection 
and transmission coefficients we used the current density, our results reduce to those obtained for 
5-type magnetic field for large values of /i, here we took \x = 1000. Moreover, the phase shift that 
characterizes the transmission amplitude is given in terms of the energy and magnetic field. 

Subsequently, we considered the case of an asymmetric double step pseudoscalar potential of the 
form W{x) = B\ [1 — Q(x + a)} + B2Q(x — a). After obtaining the general solution of the Dirac 
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equation, we matched the wavefunction at the point ±a to obtain the reflection and transmissions 
amplitudes for the negative and positive energy solutions. Different field configuration were considered 
but the interesting situation obtains when B2 = —B\ = 0.002 which enables us to recover the reflection 
and transmission coefficients for scattering of Dirac fermions by a double step pseudoscalar potential. 
Moreover, for the present potential we studied the corresponding bound states using a graphical 
solution of the transcendental equations. We reached the conclusion that there are at most two bound 
states for even parity and two for old parity solutions. 

Finally a double Gaussian type magnetic field was considered in appendix A where a second order 
differential equation for spinor components was solved. The corresponding general solution was given 
in terms of the HeunG functions. The asymptotic behavior as x — > ±00 of these wavefunctions were 
studied and enable us to extract the corresponding reflection and transmission coefficients. 
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Appendix A: Dirac fermions in a double Gaussian type magnetic field 



We consider a magnetic field normal to the graphene sheet and of a Gaussian shape. Specifically, we 
take B(x) = z \B\ sech 2 \p,(x + a)] + B2 sech 2 [/j,(x — a)]] (Fig. 12), with B\ and B2 are the strengths 
of the field. The corresponding vector potential (Fig. 12) is 



A(x) = y 



B\ B2 

— [tanh [fi(x + a)] — 1] H [tanh \p(x — a)] + 1] 

A* A 4 



(Al) 



We note that the double delta magnetic field considered in section 5 is obtained in the limit fj, 3> 1 in 
the present more general model. 

By using the same change of variable as in section 4 for ±x > 0, equation ([7]) can write in the 
following form 



4/A- 



d 



d 



dz \ dz 
for x < with q = e 2a . 



ABiqz 4B 2 qz 
{qz-lf (z-qf 



a* (qz - 1) ^i z ~q) 



0. (A2) 
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Fig. 12: The pseudoscalar potential W(x) and the associated double Gaussian type magnetic field B{x) 



(Al) for \i = 1.5 and a = 3. 



The general solution of ( A2 ) can be expressed in terms of the HeunG function [22] as follows 
^+ = Az l +{q - z) n +(-l + qz) m + HeunG (p, b+, a+, f3+,j+, 8+; ~j 



(A3) 



where we have set 



Bz- l + (q - z) n + (-1 + qz) m + HeunG (p, b 1+ , a l+ , 0+, 2 - 7+ , S+; - 



l± = (7± " 1)A n± = -«5±/2, m± = I 1 + ^ \ /2, p = <T 2 
a± = (5 2± + S 3 ±) fi~ 2 /2, ai± = (-4B 2(1) + 2 M 2 - S 2± + 5 3 ±) ^ 2 /2, 



7± = 1 + V~ 2 JfJ- 2 k 2 - /j, 2 e 2 T±B 1{2) k(i + 45 



1(2) , 



<5± = -2fi 2 B 2 (i) 



P± = [S2±S 3 ± + M 4 + (-e 2 + k 2 - A B 1{2) ) ^ 2 ± 4 B 2(1) kfi ± 45 2 2 T 45 



21 M" 



25 3± 



6± = Tl/2 



(=FmV ± 2S 2(1) ) (1 - 7 ± V 2 T mV T 2B m q 2 ^ + fc (25 1(2) g 2 + 2B 2{1) ) M 



T2 J B 1(2) 2 g 2 T 45i5 2 ] /i 5 3 ± T (Si± + (2B 1(2) + /x 2 ) 2 /x) /x 2 7± g 2 



(i 5 q 2 



h± = Tl/2 



(± M V T 2S 2(1) ) (1 - 7 ± V 2 T mV T 2B l[2) q 2 n 2 + fc (25 1(2)(/ 2 + 2B 2(1) ) fx 



T2B 1[2) 2 q 2 T mB 2 ] S 3 ± - (2 - 7± ) (S 1± + {2B l{2) + y?f /i) fi 2 q 2 



-5„-2 



s 3 ± 



with 



[2B m + [(a, + e ) fj, ± 2B 2(1) ] fi 2 [(k - e) /i ± 2fl 2(1) ] 



-2J3 2 (i) + ^ 2 7+(-) 



Si± = 

#2± = 

S 3± = ^2S 1± + ^ + {-e 2 + AB 1{2) + k 2 )fi 2 ±Ak f iB 2{1) + AB 2 {1) +4B 2 

S 4± = (2B 1 + f i 2 )(fi 2 1± + 2B 1 ). 



The asymptotic behavior of (A3) when x — > — oo corresponds to z — > and gives for and (/? 

-ifc_a: 



lim ip + = Ae ik ~ x + Be" 

x— >— oo 

lim = Ap+e ik - x + Sp-e 



-ik—x 



(A4) 
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with 



k- = \u!_ — e 



P± 



0J+ ± ik- 



uj + = k + 



2Bo 



uj- = k 



1 1 



2B 1 



In the region when x > the analogous equation to ([7]) reads 



4 ( £\ T ( 4Biqz + AB2qz 

^ Z dz \ Z dz J ^ \( q - z ) 2 (l-qzf 



2B x z 2B 2 



H{q-z) fi (1 - qz) 



+ k\ + e< 



^ = 0. (A5) 



The solution of ( A5 ) can be expressed in terms of the HeunG function as follows 



(p+ = Cz l - (q - z) n ~ (-1 + qz) m - HeunG (p, 6_, a_, /9_, 7_, ^ (A6) 
+ Dz~ l ~ (q - z) n - (-1 + gz) m - HeunG h>, 6i_ , /3_, 2 - 7 _, <y_; ^ 



where the wavefunction parameters have been defined previously. The asymptotic behavior of (A6) 
when x —> +oo, which corresponds to z — >• gives for ip + and 99" asymptotic behavior of (A3) when 
x — > — 00 corresponds to z — > and gives for y? + and 93" 



lim v?+ = Ce- ifc+a ' + De ik+X 
lim V?" = Cpe~ ik + X + J Dpe ife+x 



(A7) 



with fcj 



e 2 and p 



e 



The reflection and transmission coefficients for the positive 



and negative energy solutions are given by equation (51). With the reflection and the transmission 
amplitudes are given by r 



^ and t 

p-pi 



P p_p* ■ Figure 13 show transmission for the Dirac particle 



scattered by the pseudo-scalar potential (Al), we can see from figures 13 (a) and (d) that, the EWs, 
TWs and BS constantly appear from the energy |cj_| < |e| < |e| > |w+| and |e| < respec- 



tively. In this later case the pseudo-scalar potential (Al) behaves like a potential in well. In figures 



13 (b) and (c) the pseudo-scalar potential (Al) behaves like a potential barrier. We also notice the 
vanishing of the energy range for the EWs but the TWs and BS appear for the energy |e| > and 
|e| < \uj-\, respectively. The position of energy range for EWs, TWs and BS is affected by the value 
of the parameter p. 



- 


f , 












(b) 



■ 


f 












; f 










(d) 



Fig. 13: The transmission coefficient as a function of energy, (a) 2B\ = 2B 2 = 0.5,2.0,3.0,5.0 (left 
to right curves), (b) 2Bi = 2B 2 = -0.5,-2.0,-3.0,-5.0. (c) 2Bi = -2.0 and 2B 2 = -3.0,-6.0. 
(d) 2B 1 = -0.5 and 2B 2 = 1.0,2.0,3.0,5.0. With a = 1, k = 1 and fj, = 1. Note that, for the case 
Bi = —B 2 , we have total transmission. 
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